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Abstract 

We have derived the coefficients of the highest three 1/x -enhanced small-x logarithms of all time- 
like splitting functions and the coefficient functions for the transverse fragmentation function in 
one-particle inclusive e^e^ annihilation at (in principle) all orders in massless perturbative QCD. 
For the longitudinal fragmentation function we present the respective two highest contributions. 
These results have been obtained from KLN-related decompositions of the unfactorized fragmenta- 
tion functions in dimensional regularization and their structure imposed by the mass-factorization 
theorem. The resummation is found to completely remove the huge small-.x; spikes present in the 
fixed-order results, allowing for stable results down to very small values of the momentum fraction 
and scaling variable x. Our calculations can be extended to (at least) the corresponding a"ln^"^^ji; 
contributions to the above quantities and their counterparts in deep-inelastic scattering. 



1 Introduction 



One-hadron inclusive electron-positron annihilation, — )■ y, Z — )■ /? + X where h denotes the 
observed hadron (or a sum over all charged hadron species) and X any inclusive hadronic final state, 
is an important benchmark process in perturbative QCD which has been measured accurately over 
a wide range of centre-of-mass (CM) energies ^ The results provide crucial inputs for fit 
determinations of the fragmentation distributions (or parton fragmentation functions) Dp{x,Q^), 
see Refs. [|2]-Bll, where x represents the fraction of the momentum of the final-state parton p trans- 
ferred to the outgoing hadron h and 2^ is a hard scale, for instance the squared four- momentum q 
of the timelike virtual photon or Z-boson in the above semi-inclusive annihilation (SIA) process, 
= = s. SIA data have also provided constraints on the strong coupling constant 0. 

The theoretical description of semi-inclusive e+e annihilation is analogous to that of electron- 
hadron deep-inelastic scattering (DIS), ep — )■ e + X, via the exchange of a (spacelike) virtual 
photon or Z-boson. The SIA differential cross section can be written in terms of transverse (T), 
longitudinal (L) and asymmetric (A) fragmentation functions (timelike structure functions) |[6l, 

= ^i+-^e)f;(..e^) + isi„^eF;(x.e^) + !coseF,V.e^), (i.i) 

Here x = 2Eh/ ^/s < 1 and are the scaled energy of the hadron h and its angle relative to the 
electron beam, respectively, in the CM frame; and for photon exchange Oq = nc4Ka^/3s is the 
total cross section for Bhabha scattering times the number of colours n^- . Disregarding corrections 
suppressed by inverse powers of Q, the fragmentation functions are related to the fragmentation 
distributions by 

Fa\x,Q') = I ['^-^c^,p{zMQ'))d;:(-,Q'). (1.2) 

The coefficient functions c^^ p in Eq. (11.21 ) are known to order a^^ [TTl-fTOl, see also Ref. [fTTII . i.e., to 
the next-to-next-to-leading order (NNLO) for Ft and Fa and to the next-to-leading order (NLO) 
for Fi which vanishes for = 0. Here and throughout this article we identify, without loss of 
information, the MS renormalization and mass factorization scales with the physical hard scale Q^. 

The scale dependence of the (process-independent) final-state fragmentation distributions is 
analogous to that of the initial-state parton distributions and given by 

The (timelike) splitting functions Pj- can be expanded in powers of = as(Q^) / (47r), 

Pjf (x,as(e')) = a^pj^^^x) + alP^P^x) + alP^f^^x) + ... . (1.4) 

The leading-order (LO) and NLO contributions P^^)^ and /'(i)^ to Eq. (11.41) have been known for 
a long time [[T2] - fT6ll . A direct calculation of the NNLO corrections P^^)^ has not been performed 
so far. However, an indirect determination [[TTlfTTl . using non-trivial relations to the spacelike DIS 
case [18] and the supersymmetric limit [fT3l[T5l[T9] - l23l has been completed recently [24] up to a 
minor caveat, which is not relevant in the present context, concerning the quark-gluon splitting. 
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Eq. (11.41 ) and the corresponding fixed-order approximations to the coefficient functions (see 
below) are adequate except for l—x<^ 1 and x <^ I, where higher-order corrections generally 
include double logarithms which can spoil the perturbative expansions. Here we focus on the 
small-jc case, where the leading contributions to the N"LO splitting functions are of the form 



r,(n)T , \ c- (") 1 1 2w , (w) (n) c\ 

Pji W = Ojga] -In X + ... , flq = — ag (1.5) 
X Ca 



where 5,j is the Kronecker symbol, and Ca and Cp are the standard SU(N) colour factors, with 
Ca — fic = 3 and Cp = 4/3 in QCD. The coefficients a^"^ and the corresponding subleading 
contributions lead to corrections which are numerically far larger than the corresponding single- 
logarithmic enhancement of the analogous spacelike N"LO splitting functions governing the DIS 
case If25ti29ll : for n = 2 see Figs. 1 of Refs. [[13 and (241 . On the other hand, the all-order Mellin- 
space summation of the leading-logarithmic (LL) contributions (11.51) leads to [|30l 



^P^^{N,a,) = P^^{N,as) = ^ 1) |(l + ^7^1 ^ " 1 1^ + NLLterms (1.6) 

which can be expanded to all orders in x-space via the standard Mellin transform 

1 



M 



- In'^x 

X 



W i ^ ^ . 0.7, 



Eq. (|1.6i) corresponds to a small and oscillating function in jc-space, suggesting that the small-x 
enhancement of P^^ (^, OCs) and (x, ttg) - which is negative in the former and positive in 
the latter case, see below - is unphysical and can be removed by extending Eq. (11.61) to the next- 
to-leading logarithmic (NLL) and next-to-next-to-leading logarithmic (NNLL) small-.x; accuracy. 
Even the former extension has not been performed in the MS scheme so far, as the results of 
Ref. [31] are not given in this scheme (and consequently do not agree with the NNLO next-to- 
leading logarithms of Refs. |[r7[[24l ). For a detailed discussion see Ref. [32] where also the LL 
result for the MS transverse coefficient function Cj g corresponding to Eq. (|1.6I) has been derived. 

In this article we employ constraints provided by the structure of the unfactorized fragmenta- 
tion functions in dimensional regularization (33) and the all-order mass -factorization formula to 
derive the coefficients of the respective highest three non- vanishing logarithms for all four timelike 
splitting functions Pj^{x, tts), /, j = q, g, as well as the corresponding coefficients for both coeffi- 
cient functions for Fj, to all (in practice sixteen) orders in ttg. The derivation of the second/third 
logarithms is made possible by the NLO/NNLO fixed-order results; consequently only the highest 
two logarithms can be resummed for the longitudinal fragmentation function Fp. 

The remainder of this article is organized as follows: In Section 2 we describe the theoretical 
framework used to perform the resummation and comment on the calculations which were carried 
out using the latest version of Form and TForm f34l[35l . The resummed splitting functions 
are written down and discussed in Section 3. The corresponding results for the transverse and 
longitudinal coefficient functions are presented in Sections 4 and 5, respectively. Our findings are 
summarized in Section 6, which also provides a brief outlook to future applications and extensions. 
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2 Method and calculation 



The main quantities in our resummation are the unfactorized flavour-singlet partonic fragmentation 
functions in Z) = 4 — 2e dimensions 



(2.1) 



where the summation over z = q, g and the MS removal of {AnY and factors [36] are understood. 
The D-dimensional coefficient functions Ca,, include all non-negative powers of £ in Eq. (|2.1I) . 



(2.2) 



k=0 



Besides the fragmentation functions Ft and Ft of Eq. (|l.ll) - is a non-singlet quantity without 
l/x terms - we consider SIA with an intermediate scalar (|) coupling directly only to gluons via an 
additional term (|) G'"^G^v in the Lagrangian, where G*"^ represents the gluon field strength tensor. 
Such an interaction, suggested as a QCD trick in Ref. [37], does occur in the Standard Model for 
the Higgs boson in the limit of a very heavy top quark [38]. The NLO and NNLO quark and gluon 
coefficient functions for the resulting fragmentation function F(f have been presented in Ref. [[24|. 

The final-state transition functions Zfj^ collecting all negative powers of £ are related to the 
matrix of the splitting functions in Eq. (|1.4I) by 



Jin 2^ 



(Z^) ' = Po(as) 



dZ^ 

da^ 



[z^y' 



(2.3) 



with 



pi pi 

qq gq 

pT pT 

qg gg 



and |3z)(as 



£as — |3oag — |3i a 



(2.4) 



As we shall see from the next equation, only the leading coefficient of the four-dimensional beta 
function of QCD, |3o = 11/3 C^ — 2/3 f39l where rij: stands for the number of effectively mass- 
less quark flavours, enters the resummation of the highest three small-x logarithms. 



Eq. (|2.3I) can be solved for Z order by order in ttj. Suppressing all functional dependences, as 
already done for most quantities in the previous two equations, the first four orders read 



1 +as^yo + «s^{2^(yo-Po)Yo + ^yi} 



+ 



3/ 1 



+ a\ 



l6£3 
l24£4 

1 

+ 
+ 



(Yo - Po) (Yo - 2po) Yo + ^ (Yo - 2|3o) Yi + (Yi - Pi ) 2Yo 



(Yo-Po) (Yo-2po) (Yo-3po)Yo 



+ 3i^^ 



24£3 

1 

24£2 



(Yo - 2po) (Yo - 3po) Yi + (Yo " 3po) (Yi - Pi) 2Yo + (Yi - 2pi) (Yo - Po) 3Yo 



(Yo - 3|3o) 2y2 + (Yi - 2|3i) 3yi + (Y2 - P2) 6Yo 



+ 4i^3 



+ 



(2.5) 
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The corresponding higher-order contributions have been generated in FORM to order aj^. It is 
clear from these results that the N"LO corrections, i.e., the splitting functions up to Y„ = y*^"-* 
defined analogous to Eq. (11.41) together with Po, (3„, determine the highest n+l powers of 
1/e in Eq. (|2.5I) at all orders in tts. Keeping in mind y^oc 1/(A/^— l)^"+\ one notices that Po and 
Pq enter at NLL and NNLL small-x accuracy only. Moreover Pi is suppressed by three powers in 
1/(A'^ — 1) relative to ; hence this coefficient contributes only at the level of the fourth logarithms, 
i.e., beyond our present accuracy. 

The same considerations apply to the unfactorized structure functions (12.11) . which at N"LO 
require the coefficients c^^f^ with £ + k < n for Fj and Fi^, and £ + k < n+l for Fi in Eq. (12.21) . 
For the convenience of the reader we collect the coefficient function results which form the input 
of the small-jc resummation discussed below. The expansions about N = N —I = for read 



c 



(2,0) 
T,q 



^CFHfN^^ + fCfUfN^^ -^CftifN^^ (2.6) 



and 



4^g^^ ^ -16CfN'^ -SCfN'^ - {l2-6i;2)CFN-^ , 

c!j^f ^ leOCFCAN-"^ -^CfCaN'^ - [{^ + 161;2)CfCa + SC^]N'^ . (2.7) 
The corresponding results for F(^ are given by 

^(2^0) ^ M(^^^^^-3_^c^^^^-2_^^^^^^-l (2.8) 

and 



4'f = -SCAN-'+[C4-^^2)CA-lnf]N\ 

cl'f ^ -16CaN-^ + 6^2CaN-' , 

cl'f = mC^N-'^ - [^Cl + fCAnf-fCFnf]N-' 

- [{^^ + l6i;2)cl-^-fCAnf-^CFnf]N-\ (2.9) 
The coefficient functions for Fi are, to the lesser accuracy required in the present context, 

4;f = 2CF + 0iN) , 41^1) = OiN') , 

4^*^^ = -fCFrifN-^-SCFrifN-^ (2.10) 



and 



4|g°^ = 4CfN^^ -4CfN^ : c[]^^^ ^ SCfN^^ + SCfN^\ 

^ -64CfCaN-^ + [^CfCa-16C^]N-^ . (2.11) 
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Note that our normalizations of Cj „ and g differ by a factor of 1/2 from those in Refs. [|9l[T0l. 
Eqs. (|2.6I) - (|2.1 II) . and some contributions with a higher l-\-k used to further overconstrain the 
systems of equations discussed below Eq. (12.151) . have been obtained from the full x-space ex- 
pressions in terms of harmonic polylogarithms (HPLs) as discussed in Ref. [|40l and coded in the 
Harmpol package for Form [|34l together with Eq. (fTTTl) . 

The corresponding expressions for the NLO and NNLO splitting functions can be read off 
directly from Eqs. (13) and (14) in Ref. and Eqs. (20) - (23) in Ref. [IH. For completeness 
we finally give the small-A'^ expansions of the LO splitting functions which we need to order A'^ ^ , 

= {i-^^2)CFN^0{N^) , P^^ = luf-fufN+OiN^), 

p^^^ = 4CfN-^ - 3Cf + ICpN + 0{N^) , 

P'i^ = ACaN-' - '-^CA-lnf+ {fCA-A:,2)CAN + 0{N^) . (2.12) 

An easy way to obtain the coefficients of any desired positive power of A'^ is to transform the 
functions to A'^-space harmonic sums [|4T|. multiply by a sufficiently large power of A'^^^, transform 
back to x-space and proceed as above. Routines for the Mellin transform of the HPLs and its inverse 
are also provided by the Harmpol package. 

Inserting the A^-space small-x expansions (12.61) - (12.121) into Eqs. (12.11) - (12.51) . we obtain the 
highest three (two) logarithms for the a"£^"+^, £ = 0, 1, 2 (£ = 1, 2), contributions to Fj^k and 
Fisf.k {Fi k) to all orders in tts for which the higher-order extension of Eq. (12.51) has been coded. It 
turns out that the a" contributions to g for a = T, (|) can be written as 

Fi"g)(A^,e) = V — — ^7 TFT i^ag^ + ^B^f + e'cfg"^ + ■ ■ ■ ) (2-13) 

or 

= I x-'-'('^-'^HAi'f+eBi'f+e'ci:,"^ + ...) (2.14) 

up to terms of order (A^— 1)^, i.e., non-x^^ contributions. Eqs. (12.131 ) and (12.141) and the corre- 
sponding results for Fr^q, F(|),q and F^i given below form the crucial observation of this article. 

Focusing for a moment on the leading logarithms, Eq. (12.141) decomposes f}^, which includes 
terms of the form x ^ ln"+^^x at all orders e with 6 = 0, 1 , 2, . . . , into n contributions of the 
form 

g-2«+i_^-i-fc£ _ e-2"+^x-^[\ -k£\nx+ ^{k£f\nh+ ...] (2.15) 

with k = 2,4, ...,2n. Since Fa^ only starts at order £^", the coefficients Aa£^ in Eq. (12.141) have 
to be such that the coefficients of £^, . . . , £"^^ in the square bracket in Eq. (12.151 ) cancel in the sum 
of these n contributions. Together with the three non- vanishing coefficients coefficients of £^"+^, 
i = 0, 1, 2, in Fa\g known from the above NNLO results, we thus have an overconstrained system 

(£ n) 

of n + 2 linear equations for the n coefficients g at each order n of the strong coupling. It is 
non-trivial that all these systems have solutions, e.g., there would be no solutions if the factor of 
two in front of (n — i) in Eqs. (|2.13l) and (|2.14l) was absent, or if the sign of this term was different. 
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The situation is completely analogous for the second and third logarithms. The splitting func- 
tions and coefficient functions up to NNLO lead to n + 1 equations for the coefficients B( „ , and to 
n equations for the coefficients Q „ in Eqs. (12.131) and (12.141) . Also the latter system can be over- 
constrained at all orders except for n = 3 and n = 4, from which the corresponding contributions 
to the N^LO coefficient functions c^^ = c^-^^ and splitting functions Pj^'^ are determined. 

The decomposition corresponding to Eq. (12.131) for pj^q, a = T,^, which are suppressed by 
one power of (A^^ — 1) or \nx relative to the gluonic quantities, is given by 

fW^e) = ^ E 1 J(„ _ (<f + + e^ci.'if' + . . . ) (2.16) 

for n > 1 (there are no x^^ terms at order tts in these cases, see Eqs. (12.61) . ( 12.81) and (12.121) above). 
The missing equation, due to the lack of an e^"+^ contribution, is compensated by the absence 
of an x~^^^^ term in the decomposition. Consequently also the three coefficients written out in 
Eq. (12.161) can be determined from the NNLO quantities given above. 

We have solved the systems of equations for these coefficients and their gluonic counterparts in 
Eq. (12.131) at all orders evaluated for in Eq. (12.51) . i.e., to order a^^^. Re-inserting the results into 
these equations then determines the respective highest three logarithms in F^f"^^^ for a = T,^ 
and = q, g to all orders in e, after which the mass -factorization can be performed to this order 
in ttg. It is worthwhile to recall that, since the coefficients of £^", . . . , e^^ at order a" are given 
in terms of lower-lower quantities, this process includes a very large number of automatic checks. 
Also these steps have been carried out using Form and, for the more involved last step, TFORM. 
The resulting splitting functions and coefficient functions are presented in the next two sections. 

Analogous to Eqs. (|2.13l) and (12.161) the unfactorized partonic longitudinal fragmentation func- 
tions at all orders n>2 can be decomposed as 

- ^E ^_lJ(„_Qe (<^<?'^-) 
up to terms of order (A'^— 1)*^. Due to the additional factor of £ relative to the previous cases, 

(£ n) 

the determination of the third coefficients ■ would require the presently unknown third-order 

coefficient functions. The determination of the two highest logarithms in ci"f^ is performed in the 
manner discussed in the previous paragraph, and provides additional checks of the splitting func- 
tions determined from Fj and F^. The resulting coefficient functions are presented in Section 5. 

Like their counterparts for the large-x limit in DIS in Ref. [|42]| (the publication of the corre- 
sponding analysis of SIA is in preparation ^431), see also Refs. [I44II45L the decompositions (12.131 ) - 
(12.181 ) are inspired by and related (but not identical) to the decomposition into purely real-emission 
and the various mixed real-virtual contributions. The cancellations of, e.g., the £^^"+\ . . . , £^"+^ 
terms between the n contributions to Eq. (12.131) are thus related to the KLN theorem [|46l . 
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3 Resummed timelike splitting functions 



We are now ready to present our (mostly) new all-order small-x results. With the exception of 
graphical illustrations, we will continue to work in Mellin-A^ space. Recall that the connection 
to x-space is simple except for the coefficients of (A^— 1)^ with k > in the expansion of the 
lowest order quantities about N = I which are required for the all-order mass factorization. These 
coefficients are not included in the all-order formulae below. 

In this section we present the resummed timelike splitting functions to next-to-next-to-leading 
logarithmic (NNLL) accuracy, 

oo 

The leading log (LL) and next-to-leading log (NLL) contributions for P^^ and P^^ have the form 



gg.LL 



(-8Q)"+1 in) 
2(A^-1)2«+1 gi 



and 



^gq,NLL 



(-8)"C; 
3(A^-1)2" 

3{N-\ 



;ilCi + 2Qn^)5Wj 



\2n 



(3.2) 



(3.3) 



(3.4) 



The coefficients in Eqs. (13.21) - (13.41) have been determined to order aj^ (n = 15 in Eq. (13.11) ). and 
are given in Table [H to the tenth order in tts - for the next six orders see the text below Eq. (13.131) . 



The the highest two contributions to P^„ and P^„ can be written as 



qg,NLL 



qq 

(n)T 
qq,NLL 



{N) 



3{N-iy" qi 



and 



qg,NNL 
qq,NNL 



(N) 



9(A^-l)2«-i 



9(A^_l)2«-i 



(3.5) 



(3.6) 



(3.7) 



The coefficients in Eqs. (13.51) - (13.71) are given in Table [2] to the sixteenth order in tts, for brevity 
using a numerical form for n > 12. 

The general form and generating function for these series are known at this point (to this author) 
only for Eq. (13.21) and the non-C^ terms in the square brackets in Eqs. (13.31) and (13.41) . i.e., those 
entries that do not involve factorial denominators. A^"-* are the Catalan numbers n47[|48ll . 



(„) ^ {2n)\ ^ J_ 
n\{n + l)\ n + 1 



In 
n 



(3.8) 
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n 


4(") 


r(") 


d(«) 
^gg.2 


D(n) 
^gq,l 


^gq,2 


^gq3 


A in) 


U 


1 


1 




o 

y 








1 


1 


1 


2 


9 


— 


— 


— 


2 


2 


3 


5 


29 


1 


1 


1 


3 


5 


10 


49 
T 


100 


5 


19 

T 


11 

T 


4 


14 


35 


347 
~6~ 


357 


21 


179 
~6~ 


73 


5 


42 


126 


6353 
30 


1302 


84 


3833 
30 


1207 
30 


6 


132 


462 


11839 
15 


4818 


330 


7879 
15 


2021 
15 


7 


429 


1716 


624557 

o 1 r\ 
zlU 


18018 


1287 


444377 

O 1 A 
ZlU 


96163 

O 1 A 
ZlU 


8 


1430 


6435 


316175 
28 


67925 


5005 


236095 
28 


44185 
28 


9 


4862 


24310 


54324719 
1260 


257686 


19448 


42072479 
1260 


6936481 
1260 



Table 1: The coefficients of tlie LL and NLL small-x approximations in A^-space (|3.2I) - (13.41) for 
the timelike gluon-gluon and gluon-quark splitting functions for the first ten orders in ttj. Also 
shown (last column) are the related NLL quark-parton coefficients in Eq. (I3.5I ). 



5^gji and fi^qjj are given by EH 

in) _ f2n-l\ [n] _ f2n-l\ _ (n) (n) 

and the remaining coefficient in Eq. (|3.4I) is related to these results by 

^gq!i = ll^gg!i-24f • (3-10) 
Furthermore it is interesting to note that the last entries in Eqs. (|3.3I) and (|3.4I) have a much simpler 
difference, 

and that the quark-parton coefficients in Eq. (13.51 ) are related to the above quantities by 

Hence only one more complicated series is contained in \, \ and A ■ ; and an analytic 

ggi^ g4P 

formula for any of these quantities would lead to closed expressions for all a"^^/{N— 1)^""^ 



contributions to the timelike splitting functions. The coefficients in Eq. (13.51) for n = 10, . . . , 15 are 
(10) _ 12229277 (n) _ 136789507 (12) _ 245398487 

'i' ~ 630 ' ~ 1980 ' 'i' ~ 990 ' 

^(13) _ 16139182231 ^(u) _ 6986603759 ^(15) _ 102190158383 
~ 18018 ' 1' ~ 2145 ' ~ 8580 ' 
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n 


qg 


rW2 

qg 


qg 


qq 


r(«)2 

qq 


r(«)3 

iJrtn 


2 


2 


1 


2 


1 




_ 


3 


89 

y 


17 

y 


10 


26 


2 


8 
3 




863 


49 


253 


395 


37 


107 


4 


~6~ 


y 


~6~ 




y 


~6~ 




3028 


997 


5153 


16961 


892 


913 


5 


5 


To" 


30 


30 


ly 


To" 


6 


219389 


3574 


20701 


207263 


11807 


12617 




90 


9 


30 


90 


45 


30 


7 


6069467 


197447 


57770 


412927 


349373 


64299 


630 


126 


21 


45 


315 


35 




38066605 


2216267 


9169289 


36475945 


11537219 


6518189 


8 


1008 


360 


840 


1008 


2520 


840 




35395649 


8688247 


7255001 


143112541 


46944767 


81079091 


9 


240 


360 


168 


1008 


2520 


2520 


10 


1449188057 


595482761 


537123949 


22226523 


157729997 


59261597 


7S70 


\JJ\J\J 


^1 so 


4.0 


7100 


/ISO 


1 1 


41422167289 


17097960349 


93203672711 


120436395671 


41718615557 


8228126859 


18480 


46200 


138600 


55440 


138600 


15400 


12 


8.7381596 10" 


1.4491517 lO'' 


2.6500373 10" 


8.4902823 10" 


1.2012745 lO*' 


2.1542828 10" 


13 


3.4082509 10^ 


5.6761999 10^ 


1.0438829 10^ 


3.318678410' 


4.7804744 10^ 


8.6473779 10^ 


14 


1.3302661 10^ 


2.2242725 10^ 


4.1111839 10' 


1.2976945 10^ 


1.8985567 lO' 


3.4597524 lO' 


15 


5.1960779 10^ 


8.7203656 10^ 


1.6190693 10*^ 


5.0769752 10^ 


7.5293381 lO' 


1.3808638 10^ 



Table 2: The corresponding coefficients in Eqs. (13.51) - (13.71) for the timelike quark-gluon and 
quark-quark splitting functions to the sixteenth order in the strong coupling constant. 



The corresponding expressions for Eqs. (13.21) - (13.41) can be inferred from Eqs. (13.81) - (|3.12l) . 
The results (13.81) - (13.101) lead to the closed NLL expressions 



and 



with 



Cf=Q 



— P'^iN) 



tNLL 
Cf=0 



{(1-4^)1/2-1} i (A^- 1)2 (1+„^/Ca) 
{(l-4^)-i/2+i| (^Q + in^) as 



8CAas , tts 
and = — 

(A^-l)2 471 



(3.14) 



(3.15) 
(3.16) 



The first line of Eq. (13.141 ) and the directly related LL part of Pg^^{N) agree, of course, with the 
classical resuk (fTel) of Refs. [30|. Already at order a^^ [ElIM], the NLL second line of Eq. (|3.14l) 
is not the same as the result in Ref. [f31] which does not refer to the MS scheme, see Ref. [32] . 

The expressions for the third logarithms (NNLL for Pgi and N-^LL for iqi) are far more lengthy. 
We therefore confine ourselves here to the full analytic expressions at order a^, and present the 
higher-order coefficients only in numerical form for the case of QCD, Ca = 3 and C/r = 4/3. The 

(3)T (3)T 

leading N^l behaviour of Pgg and Pgq is given by 

= -^20Cl + ^{li^Q^ + fQV-fciC,.,} 

+ ^ CACpn} - if C^nf } + . . . (3.17) 

with ^2 = ^^/6 and 

I f{ iw ,26 r\r^r (55 22>-\^2^2 193 r'2-r V. 

+ {^M- 1)5 I ' 3" + T ^A^F -\n-- ^A^F - m^A^F^f 

+ ^ CACjrif - i CaCfu] + ^ Cjnj } + .... (3.18) 

The corresponding results for P^^ and P^'^ read 

512 00 512 

512 



+ (aT^ { (- W + I C2) Cin, + (i + i ^2) CiCn, - § Cin/ 

+ ^CACFn/-27CAnj + ^CF/ij} + ... (3.19) 



and 



qq 



(A^-l)6 y ^ ' / (A^ 
512 



-f^CACpnj^f^Cjnf^ + ... , (3.20) 

where the dots indicate terms beyond the present accuracy of the expansion in powers of 1 / (A'^ — 1 ) . 
The respective NNLL and N^LL higher-order expressions are written as 

'•IJJnlC^) = JW^^. ('""CS + 96'- C|l + „/) ,3.21) 
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and 



''S'Lm = y^^y|;r2(''''-'<"l«/ + 96"-^<3"? + 96^^^^ (3.22) 

The coefficients for Eq. (13.211) and Eq. (13.221) are given in Table 3 and Table 4, respectively. 
Here the relative normalization of the coefficients of different orders in is such that the ra- 
tios will tend to one for n — oo, if the NNLL correction have the same convergence 
properties as the LL and NLL contributions in Eqs. (13.141) and (13.151) . The present calculations 
have not been carried out to an order sufficient to definitely decide whether this is indeed the case. 

The fixed-order and resummed timelike splitting functions are illustrated and compared in 
Figs. [i|-|4]at a standard reference scale, ~ A/J, for = 5 effectively massless flavours. For 

the corresponding value tts — 0. 12 of the strong coupling constant, the expansions to order a^^ are 
sufficient, and for some of the NNLL results required, for an accuracy of 0.1% or better down to 
the lowest jc- values shown, x= 10^. An extension of the maximal order to cover one more order 
of magnitude in x is definitely feasible, but does not appear to be warranted for any foreseeable 
analyses of experimental data. 

It is clear from Figs. [T] and [2] that the available fixed-order approximations to the splitting 
functions are not reliable at x < 10^^ for the gluon-parton cases - recall Eq. (11.31) and the form 
(12.41) of the timelike splitting function matrix, which is transposed relative to the spacelike case of 
the initial-state parton distributions - and x ^ 10^ for the quark-parton cases. Obviously it is also 
insufficient to only add the previously known leading-logarithmic resummation [[30ll from order 
ttg to the NNLO gluon-quark and gluon-gluon splitting functions in Fig. [TJ On the other hand, a 
near-perfect cancellation of the strong jc-dependences is exhibited by the NNLO + NNLL results 
for xPp especially in these cases. The situation is somewhat less clear-cut for the quark-parton 
splitting functions in Fig. [2] where, as already at order but unlike the gluon-parton cases, the 
effects of the second and third logarithms have the same sign. Within the present uncertainties all 
results appear to be consistent with xPj^ ^ at .x; < 10^^. 

In Figs. [3] and m the known three fixed-order approximations are compared by their resummed 
counterparts obtained by adding the 'appropriate' higher-order resummations to the respective 
fixed-order results, i.e., forming the LO + LL (for the gluon-parton cases), NLO + NLL and NNLO 
-fNNLL combinations. The differences between the two expansions at x < 10^ are striking. 
Some questions remain due to the relatively large NNLO + NNLL corrections in Fig. [3] and the 
corresponding behaviour at x < 10^^ in Fig. HI Their answer will require the calculation of the 
fourth-order (N-^LO) splitting functions (from which the N-^LL resummations for Pgq and Pgg can 
be inferred analogously to the present calculations) which, unfortunately, is not expected in the 
near future. In the meantime the NNLO + NNLL results, and their comparison with the previous 
NLO -|- NLL resummed order, should be sufficient for practical data analysis including estimates 
of the effect of the presently unknown higher orders. 
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Table 3: The numerical coefficients of the NNLL small-x approximations (13.211) in A'^- space for the 
timelike gluon-gluon and gluon-quark splitting functions in QCD to the sixteenth order in ttj. 
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Table 4: The corresponding N-^LL coefficients in (13.221) for the timelike quark-gluon and quark- 
quark splitting functions in QCD to the sixteenth order in the strong coupling constant. 
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Figure 1 : The timelike gluon-quark and gluon-gluon splitting functions at a typical value of the 
strong coupling constant tts, multiplied by x for display purposes. Shown are the NLO and NNLO 



approximations, and the consequences of adding the leading (a" In x), next-to-leading and 
next-to-next-to-leading small-x logarithms to the latter at all higher orders in ttg. 
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Figure 2: As Fig. [H but for the timelike quark-quark and quark-gluon splitting functions, where 
the highest logarithms are of the next-to-leading logarithmic form a"^^ In^^^^x. 
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Figure 3: The timelike gluon-quark and gluon-gluon splitting functions at a typical value of tts, 
multiplied by x for display purposes. The LO, NLO and NNLO fixed-order approximations are 
compared with the small-x resummed results obtained by respectively adding the LL, NLL and 
NNLL contributions at all numerically relevant higher orders in ttj. 
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Figure 4: As Fig. [3l but for the timelike quark-quark and quark-gluon splitting functions which do 
not receive leading logarithmic (LL) corrections. Hence only two resummed curves are shown. 
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4 Resummed coefficient functions for Fj 



We now turn to the coefficient functions. For brevity, we will not discuss the (|)-exchange case here 
(beyond the respective highest logarithms which are directly related to those for Fj), as it will be 
of only theoretical interest in the near future. The corresponding results are included, however, in 
the Form file of results distributed with the arXiv version of this article. 

The moments of the small-x resummed terms of the transverse coefficient functions are 

oo 

CtAN) = I a': (5,gc W^JA^) + cg^LL W + 4?nnl W + •••)• (4-1) 

77=1 

The leading and next-to-leading logarithmic contributions for „ can be written as 



(A^_l)2n 

and 



(4.3) 



The coefficients in Eqs. (13.21) and (13.41) are given in Table |5] analytically to the twelfth order in Us 
(see the Form file for the remaining four orders) and numerically for n = 13, ... 16. In this case the 
general form and the generating function is obvious only for the leading-logarithmic coefficients 
in Eq. (Q with BSOl 

4:i = ^n(«+i) (4.4) 

k=Q 



and 

Eq. (14.51) agrees with the corresponding result of Ref. ||32l for Cj g up to a factor of two arising 
from the different normalization of this coefficient function already mentioned below Eq. (12.1 II) . 

The relations for the quark coefficient functions corresponding to Eqs. (14.21) and (14.31 ) can be 
cast in the form 



and 



2d(«) l 8/- „ dW l 8/- „ n(") 



(4.7) 



with n>2. The first sixteen coefficients in Eqs. (14.61) and (14.71) can be found in Table [6l Note the 
the faster growth of these coefficients with n, as compared to the corresponding splitting function 
results in Tables [H and [2l is largely (but only only) due to the different normalization in Eqs. (14.21) 
and Eqs. (14. 6|) . which was employed to have mainly integer coefficient in Table [B 
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Table 5: The first sixteen A^- space coefficients in Eqs. (|4.2I) - (|4.3I) for the LL and NLL small-x 
approximations to the gluon coefficient function for the transverse fragmentation function. 

As for the splitting functions, the next contributions to both transverse coefficient functions 
are considerably more complex. Since the third-order SIA coefficients functions have not been 
published so far, we give the third- and fourth-order quantities analytically. The higher orders are 
presented numerically for Ca = 3 and Cf = 4/3 below. The third-order results are given by 

4\{N) = - ^^^eOC^C, + j^^{lfclC,-iCAC,nf-f (4.8) 
~ (iV~~r)4 " ^ ^^^^ ~ (f + f ^2) CaCJ + g CACpnf - f Cfh^I + ... 
4>) = -^^fCAC.., + ^^{f^CAC.n,-|C.n/^ (4.9) 
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Table 6: As Table [5l but for the NLL and NNLL quark coefficient function in Eqs. (1461) - (14771) . 
The expansions of the fourth-order transverse coefficient functions about N = I read 

4>) = (]^390CiC,-^^{5f>ciC,-fciC.»/-f!c.C.„;} 

+ (A^ {(^ - 244&) CiC - + 18fc) CiC^ + t CiC.„, 

-5MQC|«^ + |tQC,r«/ + |lc|«/| + ... (4.10) 

and 

-|fC^n/ + |Cfnj} + ... . (4.11) 



For the coefficients of tlie tliird logarithms in Table [7] we use the notation 



4';!nnlW = jjt^- (96"C<;l„-96"-'c<;>,„, + 96"-=C<;l,„/) , (4.12) 

4'UW = (X^3^(9''"-'ci:i%-96"-'4i»;+96"-^cW^ (4.13) 

These results are illustrated in Fig.[5]for the same reference point and x-range as in the previous 
section. The situation for xc^ g and xcj ^ is largely analogous to that for the corresponding splitting 
functions xP^^ and xP^^ in Figs.[3]and|4l The NLO and NNLO fixed-order approximations (the LO 
coefficient function Cj. ^ = d{l —x) is obviously not visible in this figure) are unreliable here from 
even larger x-values than above. The small-x rise of the NNLO coefficient functions is removed 
by adding the NLL and NNLL resummations from order , leaving us with functions oscillating 
about xcj ^ ~ 0. The same behaviour, if with a considerably smaller amplitude, can be established 
down to extremely small values of x for the exactly known LL gluon coefficient function (14.51) 
already determined in Ref. Il32l . Also here it would be very interesting to known one more order in 
tts and the N^LL resummation of xcj- g. The latter, however, again requires (at least in the present 
framework) the calculation of the fourth-order contribution to the splitting function Pl^ . 

It is instructive to briefly address the impact of the (scheme-independent) LL splitting functions 
(|1.6I) and (scheme-dependent) LL coefficient functions, given in MS by Eq. (14.51 ). on the scale 
dependence of the fragmentation function and its 'gluonic' counterpart F^. This is best done by 
considering the 'timelike' physical evolution kernels Kab in Mellin space. 



Ft\ ( ^TT Ki(Sf \ f Ft 



(4.14) 



dlnQ'^ \ F^ J \ K^j J \F^ 
which are given by the matrix elements of 

K = CP'^C-' + ^{as)^C-' with C= r^'i ""^'S 1 , (4.15) 

das \c^^^ 

and the splitting function matrix (|2.4I) . In terms of powers of (A^ — 1)^^, the first term could 
be different from P^ already at leading logarithmic (a" (N — 1)^^"+^) accuracy. However, the 
relations (11.61) and Eq. (14.51 ) imply 

'■fl- = f! ^p?'') = f ' ^'^^'^ ) (4.16) 

\° ''i.u. I yo I+'-ll/ 

with c^L given by the curly bracket in Eq. (14.51) . Due to Eq. (14.161) all such contributions to the 
matrix K cancel, and the factorization- scheme independent physical kernels are correctly given by 

-^TT.LL = K^l,\A^ ^ ^ ' '^T(]),LL = ^gq,LL 7 '^(j)(|),LL = ^gg,LL • (4-17) 

A Study of the physical kernels (14.141 ) beyond the leading logarithmic accuracy could be interesting, 
but is beyond the scope of the present article. 



18 



H 










Tq,2 


^(«) 


2 


-0.0488460 








-0.0411523 








3 


-0.0052813 


-0.0044582 





-0.0543741 


-0.0020576 





4 


0.0648580 


0.2598844 


0.3689986 


0.4087098 


0.0190107 


0.0329218 


5 


0.1603366 


0.7502113 


1.1226934 


1.1073150 


0.0595507 


0.0993878 


6 


0.2804175 


1.4552223 


2.2807323 


2.0085092 


0.1178470 


0.1996540 


7 


0.4247389 


2.3722220 


3.8617425 


3.1012694 


0.1933078 


0.3345043 


8 


0.5931376 


3.5017009 


5.8830239 


4.3806502 


0.2857698 


0.5048923 


9 


0.7855611 


4.8454857 


8.3606018 


5.8442114 


0.3952535 


0.7118052 


10 


1.0020223 


6.4059897 


11.309349 


7.4907736 


0.5218662 


0.9562145 


11 


1.2425739 


8.1858750 


14.743115 


9.3198698 


0.6657575 


1.2390556 


12 


1.5072942 


10.187889 


18.674852 


11.331466 


8.2709784 


1.5612246 


13 


1.7962784 


12.414785 


23.116717 


13.525807 


1.0060672 


1.9235757 


14 


2.1096323 


14.869274 


28.080167 


15.903318 


1.2028488 


2.3269235 


15 


2.4474687 


17.554012 


33.576041 


18.464554 


1.4176262 


2.7720455 


16 


2.8099054 


20.471582 


39.614620 


21.210152 


1.6505807 


3.2596851 



Table 7: The numerical coefficients of the third small-x contributions (14.121) and (|4.13l) to the 
A'^-space gluon and quark coefficient functions for the fragmentation function Fj to order ^. 
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Figure 5: The quark and gluon coefficient functions for Ft at a typical value of ttj. Shown are 
the NLO and NNLO fixed-order approximations, and the matched LL, NLL and NNLL resummed 
results obtained (beyond LL) by adding the respective small-x terms at all relevant higher orders. 
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5 Resummed coefficient functions for Fi 



Finally we briefly present the resummed results for the longitudinal fragmentation function Fi. 
Since the NNLO (third-order) coefficient functions for this observable are not yet known, only 
the respective two highest logarithms can be resummed for both the gluon and quark coefficient 
functions. The corresponding A^-space expressions can be written as 

CL,m = I < Kc^^) + C W^LL W + 4?NNL W +■■) ^ (5-1) 

n=l 

with the gluon case given by 

(N)- (5 2) 



and 



(5.3) 



As in the transverse case, the quark coefficient functions for are suppressed by one power of 
Inx or (N — l)^^, but for n> I take the otherwise analogous forms 



and 



(n) _ 



n-4 



(5.5) 



The coefficients in Eqs. (15.21) - (15.51) are given in Tables [8] and [9l as before giving the thirteenth to 
sixteenth order in a numerical form for brevity (the exact expressions can be found in the FORM 
file distributed with this article). In this case the general formula is not even known for the LL 
coefficients which, like all other 'unsolved' series above, involve unpleasantly large prime numbers 

(7) 

early in the expansion. For instance, the prime-factor decomposition of A}^ g reads 4 • 10691. 

These results are illustrated in Fig. |6] in the same manner as those for Fj in Fig. |5] above. 
While neither of the first-order (LO) coefficient functions includes any x^^Xnx terms in the present 
case, also here the (now negative) small-x spike of both second-order (NLO) coefficient functions 
is completely removed by adding the corresponding all-order resummations of the small-jc loga- 
rithms, leaving small oscillating functions with xc^ ^ ^ at x < 10^^. 

One may expect that the small-x resummation of the longitudinal fragmentation function will 
be the first to be extended to a higher accuracy as, in contrast to the timelike splitting functions 
and the transverse fragmentation function in the previous sections, 'only' a third-order calculation 
is required for deriving the NNLO + NNLL resummation. Note, however, that already the present 
results are sufficient for the corresponding resummation of the total fragmentation function, ob- 
tained by integrating Eq. (11.11) over 0, as the coefficient functions c|^"^ are suppressed by one power 
of Injc or (A'^ — 1)^^ with respect to their transverse counterparts. 
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Table 8: Upper part: the first sixteen A'^-space coefficients in Eqs. (15.21) and (|5.3I) for the LL and 
NLL small-x resummed the gluon coefficient function for the longitudinal fragmentation function. 
Lower part: the first five (NLL and NNLL) coefficients for the corresponding quark coefficient 
function defined in Eqs. (|5.4I) and (15.51) . 



21 



n 


An) 
q 






d(«) 


^LqA 


7 


401944 


8167748 


84996 


154428 


236236 




15 


15 


5 


5 


5 


8 


6784088 


499053868 


868164 


28758068 


45616904 


o 


35 


105 


7 


105 


105 


Q 


148855862 


4294474801 


6390999 


84243073 


135167864 


y 


105 


105 


7 


35 


35 


10 


3295405924 


36585726017 


142228288 


734599784 


2115778496 




315 


105 


21 


35 


63 


11 


24496904632 


103121715842 


1061967908 


286309749296 


151010702344 




315 


35 


21 


1575 


525 




2015422894136 


85829821660568 


87708776636 


9038620655308 


6030487800584 


12 








J / / J 


/ J 


13 


4.3730248 10*^ 


2.072952610" 


2.8644635 lO'^ 


1.3405311 10'° 


2.0469392 10'° 


14 


3.302470610'° 


1.7282691 10'2 


2.169954310'° 


1.142878910" 


1.708622810" 


15 


2.503674910" 


1.4363931 10'-^ 


1.649773610" 


9.703894410" 


1.4189917 10'2 


16 


1.90453981012 


1.1906480 10 '4 


1.2582746 10'2 


8.2093343 lO'^ 


1.173626210'^ 



Table 9: Continuation of the part Table[8]for the quark coefficient function for Fi to order a^^. 
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Figure 6: The quark and gluon coefficient functions for F]^ at a typical value of ql^. Shown are 
the LO and NLO fixed-order approximations, and the matched LL (for c^g) and NLL resummed 
results obtained by adding the respective small-x terms at all relevant higher orders. 
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6 Summary and Outlook 



We have derived the all-order resummation of the highest three small-x double logarithms, 

a'^x-^ln^^-^^'-^x with £ = 0,1,2, (6.1) 

for all four flavour- singlet timelike splitting functions - with io = 2 for and P^g and io = 3 for 
Pqq and - and for both singlet coefficient functions for the transverse fragmentation function Ft 
in semi-inclusive electron-positron annihilation (SIA) - with io = 2 for Cj ^ and io= I for g - 
together with the corresponding results for SIA via an intermediate scalar (|) like the Higgs boson 
in the heavy top-quark limit. For the longitudinal fragmentation function present fixed-order 
results, which serve as input quantities for the resummation, allow only the determination of the 
highest two logarithms, i.e., £ = 0, 1 in Eq. (16.11) with io = 3 for ^ and io = 2 for g. 

The coefficients of the above logarithms have been calculated explicitly to order ^ which is 
not the highest computationally feasible order, but sufficient for numerically accurate results down 
to X = 10^^, a range in x that should be more than sufficient for all foreseeable analyses of data. 
These calculations have been performed in Mellin-A^ space, using the latest versions of FORM 
and TFORM lf34l[35l at all stages. The results agree with the leading logarithmic (LL) result of 
Refs. [30J for the splitting functions P^^ and P^g, and with the only additional result so far derived 
in the MS scheme, the recent LL contributions to the coefficient function Cj, [32J. 

The resummation has been derived by decomposing the unfactorized partonic fragmentation 
functions as,e) in dimensional regularization at any order a" into n (or n — 1 in the quark 

cases) contributions of the form 

e-^"+"<^x-^-^''^{A + Be + Ce^ + ...) with yt=l,2,...,n (6.2) 

and ^0 = 1 for a = T,^ and p = g, nQ = 2 for a = T, (|) and /? = q and for a,p = L,g, and = 3 
for = L, q, with the k = I contributions missing in the quark cases. The KLN-related cancel- 
lations between the contributions in Eq. (16.21) . together with the powers of e fixed by fixed-order 
calculations [|7l- [T7ll24l| . lead to overconstrained systems of equations for the leading logarithmic, 
next-to-leading logarithmic (NLL) [and next-to-next-to-leading logarithmic (NNLL)] expansion 
parameters A, B [and C] in the decomposition (16.21) which can be solved to (in principle) any or- 
der n. Given the large number of extra constraints and checks - including the correct predictions of 
the respective highest two small-x logarithms in the third-order timelike splitting functions I117[|24ll 
and the non-trivial all-order agreement with the known LL results Il30ll32l - there is no need for 
an additional derivation of the decomposition (16.21 ) from the structure of higher-order Feynman 
diagrams and phase-space integrations. 

Whilst the setup of the resummation is elegant and simple, most of the new results are not, as 
we have not succeeded to find the general expressions and generating functions for the resulting 
series of coefficients, with the exception of the NLL corrections to the splitting functions Cp^ P^^ 
and Pgg in the limit Cf = . The results have therefore been presented via detailed A'^-space tables 
which, hopefully, will be used for finding some of the now unknown general expressions. The 
most interesting target in this respect are the non-integer coefficients in Table [U as the solution 
of any one of these three series would be sufficient to clarify the analytic structure of all NLL 
(a"x^^ ln^"^-^x) contribution to the matrix of the timelike splitting functions. 
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The small-x resummation has a striking effect on the numerical behaviour of the splitting func- 
tions and coefficient functions in the region x ^ 10^^. All fixed-order spikes for x — )► 0, which 
dwarf their single-logarithmic counterparts in the spacelike splitting functions and deep-inelastic 
scattering (DIS) [|25] - l29ll . are removed by forming the N"LO-|-N"LL combinations of fixed-order 
and higher-order resummed results, mostly leaving small and apparently oscillating functions. This 
behaviour is qualitatively similar to the LL results of Ref. fl30l,l32il which are known in a closed 
form and thus can be evaluated down to extremely small values of x. While some theoretical ques- 
tions remain that can only be clarified by future third- and fourth-order calculations, the present 
resummation should prove sufficient for analyses of SIA data in the foreseeable future. 

We have verified that the present approach can be extended to the non-x^ double logarithms in 
the (even-A'^ based) DIS structure functions F2 and Fi (recall that there are no 'genuine' x~^ double 
logarithms in DIS; those encountered in the (|)-exchange coefficient functions in Refs. [|5T1[52|| are 
artifacts of using the heavy-top approximation outside its domain of validity). These double- 
logarithmic terms form the leading small-jc contributions in the non-singlet cases, see Refs. |[53l 
for the LL resummation of the spacelike non-singlet splitting functions; they can be relevant at 
intermediate values of x also in flavour- singlet quantities, see Ref. Il54l . The corresponding NNLL 
resummations will be presented in a subsequent publication. 

One may expect that, analogous to the large-.x: cases in Refs. [|5T1 15511 . the resummation of the 
small-.x; double logarithms can be extended to (all) higher powers of the prefactor x in Eq. (16.11) for 
the quantities considered here (and their even-N spacelike counterparts) - but not for the asymmet- 
ric fragmentation function Fa which is related to the odd-A^ structure function F3 known to receive 
additional contributions with l/nc and higher group factors [|53ll56l . We have explicitly checked 
the direct generalization of our approach to the LL and NLL x" contributions in singlet SIA for 
a = 0, . . . , 6. It works, but only for a = and even values, and with the form (16.21 ) replaced by 

£-^"+^x"-'"'{A + Be + Ce^ + ...) with k=2,...,n + l (6.3) 

which, in fact, is what one may have 'naively' expected from Refs. [9] also for the x^^ terms. 
The predictions resulting from Eq. (|6.3I) should be useful in the context of future third- and fourth- 
order calculations. Conceivably also all small-x double logarithms in the timelike and spacelike 
higher-order singlet splitting functions (and the corresponding SIA and DIS coefficient functions) 
could turn out to be 'inherited' from lower-order quantities. This issue deserves further studies 
including the case of 9\C=4- Super Yang-Mills theory addressed, for example, in Ref. [|23[|57ll . 

A Form file of our results presented in Sections 3-5 can be obtained by downloading the 
source of this article from the arXiv servers or from the author upon request. 
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